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1. Introduction 

Recent measurements of meson mass and lifetime in CDF and DO [Q] experiments are 
the first steps in the experimental research of quarkonia with open flavor. The measurement results 
are in a good agreement with the theoretical predictions for the mass 0, : 

m^^^ = 6.2756±0.0029(stat.) ±0.0025(sys.) GeV, 

mgo = 6.3000 ±0.0014(stat.) ±0.0005 (sys.) GeV, 
m^h;°'' = 6.25 ±0.03 GeV; 

as well as for the decay time ^] : 

tCDF ^ o.448lJj:[j^^(stat.) ±0.032(sys.) ps, 
tDO ^ o.475t[j:{]^^(stat.) ±0.018(sys.) ps, 

T^heor^Q 4g_^Q Q5 pS. 

Unfortunately, the experimental estimations of the cross section value were not published. 
Thus, the mechanism of Be meson production can not be understood from the obtained data due 
to poor experimental statistics, as well as due to large uncertainties in the theoretical predictions. 
Only the planed experimental research at LHC, where about 10^° events with Be mesons per year 
are expected, will improve the situation. This huge amount of events will allow to obtain the 
information on the production cross section distributions, on the decay branching fractions, and 
in some cases, on the distributions of decay products. In this research we try to fill the gap in our 
theoretical understanding of Be production and show that the study of correlations in the associative 
production of Be and D meson at LHC could be an essential information source of Be production 
mechanism. Here we study the production characteristics which weakly depend on parameters: the 
cross section distribution shapes and the ratio between S* and Be yields. 



2. Fragmentation and recombination contributions into Be production 

The Be production amplitude within the discussed approach can be subdivided into two parts: 
the hard production of two heavy quark pairs calculated in the framework of perturbative QCD and 
the soft nonperturbative binding of b and c quarks into quarkonium described by nonrelativistic 
wave function. The calculations within the discussed technique are the most simple for the process 
of Be production in the e^e^ annihilation. As it was shown in [||, ^,10], the special choice of the 



gluonic field gauge allows to interpret the Be production process as the b quark production followed 
by the fragmentation of b quark into Be meson. Thus, in the e^e^ annihilation at large energies 
the consideration of leading diagrams for the Be meson production leads the well known factorized 
formula for the cross section distribution over z = 2Eb^ / -/s: 

^ = o,-,-D{z). (2.1) 
The analytical forms of fragmentation functions for S wave states are known from 
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The relative yield of B* and Be in the e^e^ annihilation obtained within pQCD calculation 
R^g+g- = o{B*) / o{Bc) ~ 1.4. Thus the naive spin counting which fairly predicts this ratio for B* 
and B (R^+g- ~ 3) can not be applied to Be and B^ production. 

At first sight it would be reasonable to assume that for the gluonic B^ production the fragmen- 
tation mechanism is also dominant at least from transverse momenta larger than B^ mass. But as it 



was shown in [ ]11[ , |12| , |10| , |13| ] the other mechanism essentially contribute to this process practically 
all over the phase space and the fragmentation approach is valid only at transverse momenta larger 
than 5-^6 masses of Be- The total gluonic cross section predicted using full set of leading order 
diagrams essentially differ from the fragmentation approach in absolute value as well as in shape 
of interaction energy dependence. 

As it is predicted within pQCD [ pi] , [T2| ], about 90 % of the B^ mesons at LHC energies will 
be produced in the gluonic fusion (Fig. |]). Therefore in our pQCD calculations we can neglect 
the other partonic subprocess. The convolution of the gluonic subprocess cross section with the 
gluonic structure functions partially hides the differences between the pQCD predictions and the 
fragmentation approach. The predicted ratio between the hadronic cross section values is about 
of 2. Obviously, such a difference is not essential for the calculations of forth order on a^. Nev- 
ertheless, the relative yield of B* and B^ does not depend on and could indicate the produc- 
tion mechanism. Even in the kinematical region where the fragmentation model could be applied 
the value of /?hadr = 0'hadr(5*-)/o'hadr(fic ) predicted within pQCD is about 2.6 instead of 1.4 ob- 
tained within the fragmentation approach. To measure this value one need to detect B* meson 
which with unit probability decays into B^ and photon. However, it is quite difficult to detect 
such a process experimentally due to the small mass difference between B* and Be mesons [Q]: 
AM = Mb* —Mb^ = 65 lb 15 MeV. In laboratory system the maximum energy of emitted photon 

is (Omax = (^7+ \/ (7^ — 1)) where 7 is B* 7 factor, and even for B* with energy ~ 30 GeV 
Oimax ~ 0.7 GeV. Thus one can conclude that there is no certainty that the method based on sepa- 
ration of B* from Be can be used to study the Be production mechanism. This is why in the next 
chapters we suggest another way to research the Be production mechanism. 

3. The cross section distribution on the invariant mass of D and Be mesons 

Within fragmentation approach the shape of the cross section distribution over the invariant 
mass of Be and c-quark is roughly determined by b-quaik virtuality (see the diagram (3) in Fig. |T]). 
This why the distribution should be relatively narrow. Our analytical calculations confirm this 



supposition [14]. 



Here we face the problem how to transform the invariant mass of Be (B*) and c-quark Mb^+c to 
invariant mass of Be (B*) and D-meson Mb.^o- Within the fragmentation mechanism it is naturally 
to assume that the D meson takes away the total momentum of c-quark, because the production of c 
quark is the last step of emission process. Such an assumption is not obvious for the recombination 
contribution. This is why two hadronization models of c quark have been chosen: 

1. D meson takes away the total momentum of c-quark: De^o = 5(1 — z); 

2. D meson takes away part of c-quark momentum according to Kartvelishvily-Petrov-Likhoded 
fragmentation function: De^o = z^'^(l —z)- 
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Figure 3: The cross section distributions over the invariant mass of Be and D mesons calculated within 
pQCD for the process pp ^ Bc+X al = 14 TeV: without cuts (a) and for \/f > 40 GeV (b). Solid 
curves: Dc^d ~ 5(1 ~ i)\ dashed curves: Dc^d — z^'^{l — z). 

The cross section distributions over Mg^^^ are shown in Fig. |3| for the different kinematical 
regions. One can see that for the total phase space the cross section distribution looks like the 
obtained within the fragmentation approach, but the cut on interaction energy essentially transforms 
the distribution shape. It become essentially wider, whereas within the fragmentation approach it 
should remain the same. Thus one can conclude that the recombination contribution is dominant. 
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Figure 4: The cross section distribution over the 
cosine of angle between the directions of motion 
of Be and D mesons predicted within pQCD for 
the process pp Be +X at = 14 TeV. 



Figure 5: The dependencies of averaged c quark 
energy on Be meson energy are represented for 
the different cuts on ^-jet transverse momenta. 



4. The angle correlations and the decay length ratio in the associative production of 
Be and D mesons. 

As it is shown in Fig. ^, the cross section distribution on cosine of the angle between the Be 
meson and D meson has the sharp maximum at cos 6 ~ 1. Moreover, about a half of Be mesons 
is associated by the D meson moving in the close direction: 6 < 26". Therefore, D meson can be 
used to detect Be meson. 

In the Fig. § the dependencies of averaged c quark energy on Be meson energy are represented 
for the different cuts on b-jet transverse momenta. It can be concluded from these plots that at any 
Z?-jet transverse momentum 

{Ee)>l.2EB,. (4.1) 

For D meson we obtain that 

(£s)> 0.7 -1.2^5^ (4.2) 

for De^D = — z) and De-^o = 5(1 — z), correspondingly. 

The decay lengths depend on particle energies and lifetimes as follows: 

ih) - ^cTo Ib, ^ ^cTd. (4.3) 



Taking into account that Tq / , ~ 2 we obtain: 

{hai-D 



> 5. (4.4) 
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This value should be compared with the fragmentation model prediction: 




(4.5) 



5. Conclusions 

The following conclusion can be drawn from the performed calculations: 

1. The cross section distribution over the invariant mass of Be and D meson depends essentially 
on kinematical cuts and can be used to research Be production mechanism at LHC. 

2. In many cases the Be and D mesons move in close directions. It could be useful to detect Be 
meson. 

3. The energies of Be and D mesons are comparable. The decay length of D meson by more 
than 5 times larger than the decay length of Be meson. The experimental research of the ratio 
between Be meson and D meson energies could shed light on Be production mechanism. 
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